The /3-relaxation dynamics of a simple liquid 



Q\ . Tobias Gleim and Walter Kob 

qv , Institute of Physics, Johannes Gutenberg- University, Staudinger Weg 7, D-55099 Mainz, 

Germany 

a 

H- 8 

O ' Abstract 
m ' 

— ■ We present a detailed analysis of the /3-relaxation dynamics of a simple glass 

qj ■ former, a Lennard-Jones system with a stochastic dynamics. By testing the 
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various predictions of mode-coupling theory, including the recently proposed 
corrections to the asymptotic scaling laws, we come to the conclusion that in 
this time regime the dynamics is described very well by this theory. 
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In the last few years it has been demonstrated that mode-coupling theory (MCT) |1| is able 
to describe many aspects of the relaxation dynamics of supercooled liquids. In particular the 
theory is able to explain on a qualitative level, and for certain systems even on a quantitative 
one, phenomena like the non-Debye behavior of the a-relaxation process, the wave-vector 
dependence of the Lamb-Mossbauer and Debye- Waller factors, and why quantities like the 
viscosity or the a-relaxation times show an anomalously strong temperature dependence of 
their activation energy in the vicinity of T c , the so-called critical temperature in the theory. 

Apart from the existence of T c , the most important predictions of MCT deal with the so- 
called /3-relaxation process which is proposed to exist in the supercooled regime on the time 
scale between the microscopic relaxation at short times and the a-relaxation at long times. 
The /3-regime is readily seen if a time correlation function </>(£), such as the intermediate 
scattering function, is plotted versus the logarithm of time. In the supercooled regime 4>(t) 
will show at intermediate times a plateau, and the relaxation dynamics of the system on 
the time scale at which <p(t) is close to this plateau is the /3-regime. The reason for the 
existence of this plateau is that on this time scale the particles are trapped in the cages 
formed by their surrounding neighbors. Hence the predictions of the theory regarding the 
/^-regime deal with the details of the dynamics of the particles in these cages. Some of these 
predictions have already been confirmed by various experiments on colloidal suspensions and 
molecular liquids 0. 

In the past it has been shown that apart from experiments also computer simulations 
are a very useful tool to probe the dynamics of supercooled liquids [|| . Because simulations 
allow the investigations of observables which in real experiments are hard to measure, as 
e.g. the dynamics at large wave-vectors or cross-correlation functions, they permit to make 
more stringent tests of theoretical concepts and thus are a valuable addition to experiments. 
Results of such tests have, e.g., been done for soft sphere systems ||, Lennard- Jones mod- 
els || , water || , and polymers [|7| . The result of these tests was that the theory is indeed 
able to give a good description of the relaxation dynamics of these systems. What these 
simulations have, however, not been able to address so far are several important predictions, 
discussed below, of MCT about the relaxation dynamics in the /3-regime. The main reason 
why these predictions have not been tested was that they are supposed to be valid only 
very close to the critical temperature T c of MCT, and that close to T c the relaxation times 
of the system are usually so large that it is very hard to equilibrate the system within the 
time span accessible to a computer simulation (but easily reachable in a real experiment). 
If the predictions of the theory are tested at slightly higher temperatures, where the system 
can be equilibrated even in a computer simulation, the strong interference of the micro- 
scopic dynamics of the system with the /3-relaxation process will spoil the analysis, because 
of the lack of separation of time scales, and stringent tests will almost be impossible. In 
order to overcome these problems we have recently investigated the relaxation dynamics of 
a Lennard- Jones system in which the particles move according to a stochastic dynamics || . 
This dynamics leads to a strong damping of the microscopic dynamics and hence it becomes 
finally possible to test the predictions of MCT about the /3-regime and in this paper we 
report the outcome of these tests. 

The model we investigate is a 80:20 mixture of Lennard- Jones particles with mass m. In 
the following we will call the two species of particles A and B. The interaction between two 
particles of type a and /3, with «,/?€ {A, B}, is given by V a/ 3(r) = 4e Q/ 3[(cr Q , j a/r) 12 — (a^/r) 6 ] 



with caa = 1-0, cr aa = 1.0, cab = 1-5, <tab = 0.8, 6bb = 0.5, and <tbb = 0.88, and a cut-off 
radius of 2.ba a p. In the following we will always use reduced units with <taa and 6aa the unit 
of length and energy, respectively (setting the Boltzmann constant ks equal to 1.0). Time 



is measured in units of J <7 AA m/48eAA- The volume of the simulation box is kept constant 
with a box length of 9.4. The dynamics of the system is given by the stochastic equations 
of motion 



rnr 



* + V i XXa(ki-r;l) = -C*j + »&(*) • (!) 
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Here rjj (t) is a gaussian distributed white noise force with zero mean. Because of the fluc- 
tuation dissipation theorem, the magnitude of rjj(t) is related to £ by (i]j(t) ■ r]i(t')) = 
6kBT(5(t — t')5ji. We have used a value of ( = 10, which is so large that the presented 
results for the dynamics do not depend on ( anymore (apart from a trivial change of the 
time scale). Equations ([!]) were solved with a Heun algorithm with a time step of 0.008. 
The temperatures investigated were 5.0, 4.0, 3.0, 2.0, 1.0, 0.8, 0.6, 0.55, 0.5, 0.475, 0.466, 
0.452, and 0.446. At the lowest temperature the length of the run was 4 x 10 7 time steps. 
This length is not sufficiently long to equilibrate the sample. Therefore we equilibrated the 
system by means of a Newtonian dynamics for which we have found that the equilibration 
times are significantly shorter ||. Thus all the correlation functions shown in the present 
work are equilibrium curves, even if they do not decay to zero at long times. In order to 
improve the statistics of the results we averaged at each temperature over eight independent 
runs. 

In the following we will review some of the predictions of MCT about the dynamics in 
the /3-regime and will denote by <f>i(t) an arbitrary time correlation function which couples 
to density fluctuations. (Here the index I is just used to distinguish between different 
correlators.) As stated here, the predictions are valid only for temperatures slightly above 
the critical temperature T c of MCT. More general results can be found in Refs. |1]||. 

MCT predicts that in the /^-region any correlation function 4>i(t) can be written as 

Mt) = ft + hic a g4t/U) , (2) 

where the temperature independent constants ff and hi are called critical nonergodicity 
parameter and critical amplitude, respectively. The quantity c a is given by 



with a = C(T c -T) , (3) 

where C is a constant. The function g_ is independent of (pi and depends only on the so- 
called "exponent parameter" A which can be calculated from the structure factor M. This 



calculation has been done for the present system and a value of A = 0.708 was found [jlO 
Hence in our case A is not a fit parameter. Once A is known, the function g_ can be calculated 
numerically. 

The quantity t a in Eq. (^) is the time scale of the /9-relaxation and is given by 

t„ = VM 1/aa , ( 4 ) 

where to is a system universal constant, and the exponent a can be calculated from A and is 
in our case a = 0.324 flQfl . Hence, according to MCT, in Eq. (Q) only the time scale t a and 
the prefactor hic a depend on temperature. 



MCT also predicts that Ti(T), the time scale of the a-relaxation, depends on temperature 
like 

n = Tir, T = t /\a\\ with 7 = l/2a + l/26 (5) 

where Ti is independent of temperature and the exponent b can also be calculated from A |I[ 
and is for our system b = 0.627 [[HJ. Thus we have 7 = 2.34. 

Having presented some of the predictions of MCT about the /3-relaxation we can now 
check how well they agree with reality. For this we calculated from the simulation F a ^(q, t) 
and F*(q,t), the coherent and incoherent scattering functions for wave- vector q, respec- 
tively. These time correlation functions were then fitted in the /3-relaxation regime with the 
functional form given by Eq. (§), where ff, hic a , and t a were fit parameters. This fit was 
first done for the lowest temperature (T = 0.446). For the fits at the higher temperatures 
the value of ff was kept fixed to the one of T = 0.446 in order to avoid that the fits give 
some effective time scales t a and prefactors hic ff . In the inset of Fig. [J] we show the results 
of such fits and it can be seen that the range over which the /3-correlator describes the data 
increases with decreasing temperature, as predicted by MCT. 

From Eqs. @ and (Q) it follows that according to MCT a plot of t~ 2a versus temperature 
should give a straight line and that this line should be independent of the correlator <fti . In 
Fig. |l| we show such a plot where we have used for <f>i the functions F*(q,t) for the A and 
B particles and the function F AA (q,t), for two wave-vectors: q = 7.20 and q = 9.61, which 
correspond to the location of the maximum and first minimum in the structure factor for the 
A-A correlation || . From this plot we see that at low temperatures the different curves are 
indeed close to straight lines and collapse quite well onto a master curve. Hence we conclude 
that these two predictions of MCT work well for our system. Also included in the figure is 
a linear fit to the data for F A (q, t) for q = 7.2. This fit intercepts the temperature axis at 
T m 0.432, which according to MCT should be T c . This estimate of the critical temperature 
is in excellent agreement with the one of Ref. ||, where T c = 0.435 was found. 

We also mention that we have found that the square of the prefactor in Eq. (|H), hic a , 
shows a linear dependence on T, (hic^) 2 = \a\, and vanishes at T c , which follows from 
Eqs. (0) and (§) and is hence in agreement with MCT. The test of Eq. (^ is equivalent to 
the test of the relation between hiC a and t a , which according to the theory, Eqs. ©,@ and 
(|), should be 

hic a (xt- a . (6) 

In Fig. [| we plot h\c a versus l/t a in a double logarithmic plot for the same correlators 
discussed in Fig. [I]. We see that the different curves can be approximated reasonably well 
by straight lines with a slope a (bold solid line in the figure). Therefore we conclude that 
also this prediction of the theory seems to work satisfactorily well. 

From Eqs. (|J) and Eqs. (J5|) it follows that also the a-relaxation time 77 should show a 
power-law dependence on t CT , i.e. 

riccr,^ . (7) 

Thus this equation expresses the surprising prediction of MCT that two diverging time 
scale exist in supercooled liquids, namely 77 and t a . Whether this is indeed the case is tested 



in Fig. |3] where we plot rf 1 versus t^ 1 for the usual correlators in a double logarithmic plot. 
We see that the different curves are indeed close to straight lines and that the slope is very 
close to the theoretical value, bold straight line. Thus we confirm also this prediction of the 
theory. 

As already mentioned above in the context of Eq. (Q), according to MCT the whole time 
dependence of <f>i(t) is given by the /—independent function g_(t/t a ). In order to test this 
prediction we can introduce a function Ri(t) as follows: 

RM - Mt") - W) ' (8) 

Here t' and t" are arbitrary times in the j3- relaxation regime (if 7^ t"). From Eq. (0) it 
follows immediately that in the /3-regime the function Ri(t) is independent of the correlator, 
i.e. of I. To see whether this is indeed the case we have considered the correlation function 
discussed in the context of Fig. [I] and in addition the coherent and incoherent scattering 
function for several other wave vectors and also the cross correlation function F AB (q,t) at 
different q. This gave us a total of 36 correlation functions which are shown in the upper 
inset of Fig. [| (at T = 0.446). For each of these functions we determined the corresponding 
Ri{t), choosing for tf and t" a value around 200 and 15000, respectively. In the main figure 
of Fig. |] we show the different Ri(t) and we see that in the /3-regime they do indeed collapse 
onto a master curve. That such a collapse is not a trivial result can be concluded from the 
observation that outside the /3-regime the different curves show a strong dependence on I, 
at short as well as at long times. 

Equation (|2|) is the prediction of the theory about the leading asymptotic behavior for 
the time and temperature dependence of a generic correlator. Very recently the next order 
corrections to this behavior have been calculated |IT| and these corrections can now be used 



to do more checks on the validity of the theory. In Ref . [I l| it has been shown that in the 



early /3-relaxation regime, i.e. for to "C t <C t a , the correlator can be written as 

Mt) = fi+hi(to/t) a {l + [K l + A}(t /tr} . (9) 

Here A is a /—independent constant and the constant K~i depends on I but not on temper- 
ature. In the late /3-regime, for which t a -C t <C Tj, the correlation function is predicted to 
behave like 

Ut) = fi-hi{t/r) b {l-K l {t/r) h } . (10) 



The mentioned corrections are the second terms in the curly brackets in Eqs. (f|) and (|10|). 



The important result about these equations is that the / dependent part of the correction, 
i.e. Ki, is the same. Using this fact it is simple to show the following: Calculate the ratio 
Ri(t) from Eq. (|8|) for various correlators and plot these Ri(t) versus the logarithm of t. 
Draw two vertical lines at times that are a bit shorter and a bit longer than the times where 
the asymptotic expression, Eq. (Q), holds. Start to label the correlators from top to bottom 
in the order they intersect the vertical line at short times and call this number i. Determine 
the position j at which curve i intersects the vertical line at large times, where the counting 



is again done from top to bottom. Thus this gives a function j(i). From Eqs. @ and (10) 



it then follows that j — i. Or to put this in other words: the first (second, . . . ) curve 



that intersects the left vertical line is also the first (second, . . . ) curve to intersect the right 
vertical line. 

We have done the described procedure by using vertical lines at t = 3 and t = 10 5 (bold 
vertical lines in Fig. |4j). The function j(i) we find is shown in the lower inset of Fig. |j. 
We see that, despite the scattering present in the data, a clear increasing trend which is 
compatible with a straight line with unit slope can be seen, thus giving also support for the 
validity of this prediction of the theory. 

We thus can conclude from the present work that many of the predictions that mode- 
coupling theory makes for the /9-relaxation can also be tested in computer simulations. As we 
have shown in this paper the outcome of such tests for the Lennard- Jones system considered 
here is that the theory is indeed able to give a self consistent picture of the dynamics of this 
simple glass-former in the /3-relaxation regime. 
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FIG. 1. Main figure: Check of the validity of Eq. (Q) for various correlators (see labels of 
curves). MCT predicts a straight line which intercept the T— axis at T c = 0.435. The bold straight 
line is a linear fit to the open circles. Inset: Time dependence of F s (q,t) for the A particles at 
q = 7.2 (symbols) and the fitted /3-correlators (solid curves) for all T < 0.8. 
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2. Check of the validity of Eq. (0) for various correlators. The prediction of MCT is a 
line with slope a (bold straight line). 
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3. Check of the validity of Eq. (0) for various correlators. The prediction of MCT is a 
line with slope 1 + a/b (bold straight line). 
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FIG. 4. Main figure: Time dependence of the ratio given in Eq. (||), demonstrating the validity 
of the factorization property [ Eq. (g)]. T = 0.446. The correlation functions 4>i(t) are shown in 
the upper right inset. See text for the discussion of the lower left inset and the two vertical lines 
in the main figure. 



